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ABSTRACT 

Context. Galaxy formation in the current cosmological paradigm is a very complex process in which inflows, outflows, 
interactions and mergers are common events. These processes can redistribute the angular momentum content of 
baryons. Recent observational results suggest that disc formed conserving angular momentum while elliptical galaxies, 
albeit losing angular momentum, determine a correlation between the specific angular momentum of the galaxy and the 
stellar mass. These observations provide stringent constraints for galaxy formation models in a hierarchical clustering 
scenario. 

Aims. We aim to analyse the specific angular momentum content of the disc and bulge components as a function of 
virial mass, stellar mass and redshift. We also estimate the size of the simulated galaxies and confront them with 
observations. 

Methods. We use cosmological hydrodynamical simulations that include an effective, physically-motivated Supernova 
feedback which is able to regulate the star formation in haloes of different masses. We analyse the morphology and 
formation history of a sample of galaxies in a cosmological simulation by performing a bulge-disc decomposition of the 
analysed systems and their progenitors. We estimate the angular momentum content of the stellar and gaseous discs, 
stellar bulges and total baryons. 

Results. In agreement with recent observational findings, our simulated galaxies have disc and spheroid components 
whose specific angular momentum contents determine correlations with the stellar and dark matter masses with the 
same slope, although the spheroidal components are off-set by a fixed fraction. The average angular momentum efficiency 
for the simulated discs is 77 ~ 1 while bulges is 77 ~ 0.10 — 0.20. For the simulated sample, the correlations found for 
the specific angular momentum content as a function of virial mass or stellar mass are found not to evolve significantly 
with redshift (up to a ~ 2). Both dynamical components seem to move along the correlations as they evolve. The total 
specific angular momentum of galaxies occupy different positions filling the gap between pure rotational-dominated and 
dispersion-dominated systems. The scaling relations derived from the simulated galaxies determine a similar relation 
with the virial radius, in agreement with recent observations. 

Conclusions. 

Key words, galaxies: formation, galaxies: evolution, cosmology: dark matter 


1. Introduction 

The understanding of the morphology of galaxies and its 
relation with their physical properties has evolved during 
the last decades as new insights are provided by observa¬ 
tions. It is well-known that the angular momentum con¬ 
tent correlates wi th the mor phology of galaxies along the 
Hubble sequence. iFalll (|l983ll has found evidences suggest¬ 
ing that the specific angular momentum content of both 
disc and elliptical galaxies holds a clear correlation with 
the stellar mass, but wi th the elliptical systems it shifts 
to lower values. Recentlv iRomanowskv fc Falll (|2012ll have 
analysed a larger database of galax ies w i th dif ferent mor¬ 
phologies, revisiting and confirming iFalll (|l983ri ’s findings. 
These authors estimated the specific angular momentum of 
the stellar component js using kinetic and photometric data 
and the stellar mass obtained from near infrared luminosi¬ 
ties. Using this larger database, they found that not only 
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disc-dominated galaxies follow this relation but also ellipti¬ 
cal ones. However, the latter has lower js at a given Mstar- 
The authors suggested that observed galaxies of all mor¬ 
phological types lie along almost parallel sequences with 
exponents a ~ 0.6 in the js — Mstai diagram, but with an 
off-set between late and early type. More recently, citet- 
Fall2013 showed that galaxies position in the js — Mstar 
diagram is correlated with their position in the disc-to- 
bulge-ratio-Mstar plane. They concluded that the descrip¬ 
tion of galaxies in terms of js — Mstar may be regarded as 
a physically motivated alternative to the traditional mor¬ 
phological classification. The detected js — Mgtar relation 
for spheroidal galaxies suggests an an gular momentum effi¬ 
cienc y of 10 — 20% regardless of Mstar ( Fall fc Romanowskvl 

IMl- 

In a hierarchical clustering universe, protogalaxies ac¬ 
quire their angular momenta via t idal torques due to small 
density fluctuations (|Peebleslll969h . As baryons are initially 
drawn by the same tidal torques as the dark matter haloes. 
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they are expected to acquire similar specific angular mo¬ 
menta. The specific angular momentum of the dark mat¬ 
ter haloes is supposed to grow as . Hence, if baryons 
collapse conserving their angular momentum, they should 
also determine a s i milar relation. Based on this hypothesis, 
iFall fc Efstathioul (|198Clll postulated the so-called standard 
disc model (SDM). This model was able to explain the ob¬ 
served correlation betw een the angular momentum of discs 
and the stellar masses (|Mo et al.lll998ll . 


Several numerical works have shown that if a disc com¬ 
ponent is able to assemble then it does so conserving 
its specific angular momentum. Although within the hi¬ 
erarchical scenario for galaxy formation, this might im¬ 
ply a combination of processes which could contribute 
to retain or balance the angular momentum content so 
that discs are formed as expected from the SDM, (e.g 


that discs are lormed as expected trom the bUM, (e.g . 
Agertz et al. l20Hl iBrook e t al.ll20 1ll: De Rossi et al.ll2012t 


Padilla et al.l 12014 iNelson et ai.l 1201511 . In particular. SN 


feedback plays a critical role in determining the fraction of 


Dommffuez-Tenreiro et al. 1998 

; Scannaoieco et al. 

2009 

Guedes et al.1120111: ISales et al.l 

20121: IPedrosa et all 

2014 

Laeos et al.l 1201511. Elliotical galaxies are expected to be 


affected by major mergers or repetitive minor ones. This 
tends to decrease the angular momentum content of the 


sion( 

e.g. Meza et al.ll2003l: iParrv et al.l 

20091: iHookins et al. 

201C 

;lBournaud et al.ll2011l: Sales et al. 

I 2 OIOI: IPedrosa et al. 

2014 

). Internal disc instabilities could also contribute to 


shape bulges and elliptical galaxies (e.g. lEfstathiou et al.l 
1198211 . 


As the final morphology of a galaxy is directly linked 
to its angular momentum content, there must also be a 
connection with its characteristic size: when angular mo¬ 
mentum is cons erved, more exte nded structures are able to 
form. Recently, lKravtso\^ (l201.8ll has found that character¬ 
istic sizes of stellar and gas distributions in galaxies scale 
approximately linearly with the virial radius derived from 
the abundance matching ansatz. This author found that the 
rel ation is in agree ment with expectations from the model 
of iMo etldl (1199811 . Interestingly, he also found that this 
relation is present not only for late-type galaxies but also 
for early-type systems. From these findings, he concluded 
that the angular momentum content plays a crucial role in 
determining the sizes of galaxies of all morphological types, 
not only of those dominated by rotation. 


In a recent paper, we have made a thorough study of the 
relation betwee n the morphology of galaxies and their for¬ 
mation history. IPedrosa et al.l (|2ni4ll . The key role played 
by the angular momentum content of galaxies was studied. 
We found that when a disc component is able to form it 
does so following the expected relation from the standard 
disc for mation model in agreement with previ ous works 
(see also lGovernato et 11111201(11 : ISales et al.ll2012ll . More in¬ 
terestingly, in this simulation, spheroidal components are 
shown to follow a similar correlation but off-set with re¬ 
spect to that of rotational supported systems. It became 
clear from this work that SN feedback played a major role. 


In this paper, we analyse the angular momentum con¬ 
tent of the disc and bulge com ponents of the b est re solved 
galaxies in a A-CDM scenario. iPedrosa et al.l (120141 1 anal¬ 
ysed a sister simulation (S230A) from the Fenix Project 
with a stronger feedback parameter. We explore a run 


(S230D) with a slightly weaker SN feedback and a differ¬ 
ent chemical distribution of the synthetized metals, which 
better reproduces observational results of the metallicity 
gradients (Tissera et al. 2015, submitted). This is a very im¬ 
portant point which hel ps to set constraints on the strength 
of the SN feedback fe.g. lGibson et al.ll2013ll . The simulated 
galaxies have also larger D/T stellar ratios compared to 
those measured in S230A. We also follow them back in time 
to establish the possible existence of a redshift dependence. 

This work is organized as follows: In Section [21 we de¬ 
scribe the initial conditions and free parameters used to run 
the analysed simulation and the methods used to identify 
and decompose galaxies. In Section|21 we analyse the angu¬ 
lar momentum content of our galaxies and their dynamical 
components. In Section |31 we study the characteristic sizes 
of the galaxies, and the connection with their angular mo¬ 
mentum content. Finally, in Section |5l we summarize our 
main findings. 

2. Numerical Experiments 

We used a ve rsion of the code GADGET - 3, an update o f 
GADGET-2 ([SDringel fc Hernouistl 120031 ISpringell l2005ll . 
optimized for massive parallel simulations of highly in- 
homogeneus systems. This version of GADGET-3 includes 
treatments for metal-dependent radiative cooling, stochas¬ 
tic star formation (SF), chemical enrichment, and the 
multiphase model for the interstellar medium (ISM) and 
the Superno va (SN) feedback scheme of IScannanieco et'HI 
(120051 1200^ . The SN feedback model is able to success¬ 
fully trigger galactic mass-loaded winds without introduc¬ 
ing mass-scale parameters. 

Our code considers energy feedback by Type II (SNII) 
and Type la (SNIa) Supernovae. These processes are 
grafted int o a multiphase m o del fo r the ISM, both de¬ 
veloped bv IScannanieco et al.1 (l2006ll . The energy released 
by each SN event is distributed between the cold and 
the hot phases of the simulated ISM. The adopted mul¬ 
tiphase model allows the coexistence of diffuse and dense 
gas phases. In this model, each gas particle defines the cold 
and hot phases by using a local entropy criteria. This def¬ 
inition allows particles to decouple hydrodynamically from 
particular low-entropy ones if they are not part of a shock 
front. Each cold gas particle has a reservoir where injected 
SN energy is stored until they fullfil the conditions to join 
their local hot phase. It is worth mentioning that this SN 
feedback scheme does not include parameters that depend 
on the global properties of the given galaxy (e.g. the total 
mass, size), thus making it suitable for cosmological simu¬ 
lations where systems with different masses have formed 
in a complex way. Our code also includes t he ch emical 
evolution model developed bv iMosconi et al.l (|200ll1 . SNII 
and SNIa contribute with chem ical elements estimated by 
ado pting the chemical yie lds of IWooslev fc Weaver] (1199511 
and llwamoto et ahi (|l999ll . respectively. 

As mentioned before, for this study, we used one of the 
Fenix Project’s simulations (S230D). These share the same 
initial conditions but assume different parameters for the 
SF and SN feedback sub-grid physics. The initial conditions 
are consistent with the concordance model with = 0.7, 
= 0.3, Oh = 0.04, a normalization of the power spec¬ 
trum of as = 0.9 and Hq = lOOh km s“^ Mpc~^, with 
h = 0.7. The simulated volume represents a box of 10 
Mpc h“^ comovil with 2 x 230^ initial particles, achieving 


2 

































































































Pedrosa et al.: Angular momentum evolution for galaxies 


a mass resolution of 5.9 x Mq and 9.1 x 10®h“^ Mq 

for the dark matter and initial gas particles, respectively. 
The gravitational softening is 0.5 kpc h“^. We acknowl¬ 
edge the fact that the initial condition represents a small 
volume of the Universe. It was chosen to represent a typ¬ 
ical region of the Univ erse with no large stru cture nearby. 
Using the S230A run, iDe Rossi et al.l ( 2013li verified that 
haloes’ growth is well-reproduced by confronting the re- 
sul ts with those obtaine d from the Millennium simulation 
by Fakhouri et al. ( 201011 . The S230A was also analysed by 
iPedrosa et al.l ( ^14ll to study the morphology and angu- 
lar momentum con tent of the simulated galaxies. Recently, 
lArtale et al.l (l2015ll introduced the effects of high-mass X- 
ray binary feedback in the very early Universe in a coun¬ 
terpart run of S230D. 

Simulation S230D used a higher gas density thresh¬ 
old for star formation and a lower energy per SN event 
(0.7 X 10^^ erg) tha n in S230A (e.g. iDe Rossi et al.l [201^ 
IPedrosa et al.ll20T^ . Simulation 230D injects a larger frac¬ 
tion of new synthetized chemical elements into the cold 
phase, while the fraction of SN energy is equally distributed 
between the cold and hot phases (S230A injected equal frac¬ 
tion of energy and chemical elements into the cold and hot 
phases close to the SF region). This different combination 
of SF, SN feedback parameters and metal distribution pro¬ 
duced a good representation of chemical abundances on the 
disc components (Tissera et al., submitted). 


2.1. The simulated galaxy sample 

Our galaxy sample was built up by using a Friends- 
of-Friends algorithm to identify the virialized structures. 
Then, the substru ctures were discern ed by using the 
SUBFIND scheme (|SDringel et al.l[200Ill . which iteratively 
determines the self-bound substructures within the virial 
radius, r 2 oo- For illustration purposes, an example of 
a late type galsLxy of our simulated catalog is shown 
in Fig. [T1 (upper panels). 

To diminish resolution issues, we restricted our study 
to systems with a total number of particles greater than 
2000 in the main galaxy (which corresponds to Mvir 
lO^^'^M©). We also tested that the estimated relations are 
conserved when this particle number limit is increased to 
- 10000 . 

In order to classify our simulated galaxies into early- 
type and late-type systems, we estimated the parameter 
e of the star (gas) particles defined as e = Jz/</z,max(F^), 
where is the angular momentum component in the di¬ 
rection of the total angular momentum, and Jz,ma,x{E) is 
the maximum over all particles at a given binding en¬ 
ergy, E. A particle on a perfect prograde circular orbit in 
the disc plane should have e = 1. We considered parti¬ 
cles with e higher than 0.5 to be part of a disc. The par¬ 
ticles that did not satisfy this requirement were considered 
to belong to the spheroidal component. By inspecting the 
e — E plane, we checked whether the adopted limits were 
suitable to individualize the rotationally supported compo¬ 
nents. This t hreshold is less rest rictive than those usually 
adopted (e.g. iTissera et al.ll201^ , so if there were stars in a 
thick disc, they would contribute by decreasing the signals 
for angular momentum conservation, for example. We also 
distinguished between the central spheroid and the stel¬ 
lar halo components d e pendi ng on their bounding energy 
following iTissera et al.l (1201211 . In Fig. [T] (lower panels) we 


display the disk and central spheroid components for the 
galaxy previously shown, after the descomposition. 

Therefore, for each simulated galaxy, we determine the 
disc (i.e.rotational-dominated) and the central spheroids 
components (i.e.dispersion-dominated) based on these dy¬ 
namical criteria. This procedure was also applied to galaxies 
identified at z = 1 and z = 2 in order to study the evolu¬ 
tion of the angular momentum content in the galaxies, their 
central spheroids and disc components and the dark matter 
haloes. 


3. Angular Momentum 

To calculate the properties of the simulated galaxies, we 
adopted the half-mass stellar ratio (rnivi)) defined as the 
radius that enclosed 50% of the total stellar mass. This 
characteristic radius was chosen in order to make an ade¬ 
quate comparison with observational results. The half-mass 
radius is usually adopted to estimat e gala xy properties in 
observations. In the case of iKravtsovI (|2013f) for example, he 
uses the stellar half mass radius to estimate his Rhm — Rvir 
relation. We estimated the ratio between the angular mo¬ 
mentum of the stars (gas) in the discs (Jd) within thm 
and the corresponding of the dark matter haloes within the 
virial radius (Jh), jd = </d/Jh, as a function of the ratio 
between the corresponding ma sses, mn = Mn/M rE]. We fol¬ 
lowed the classical notation of iMo et (Il998ll to define jd 
and md. Similarly, we defined js and mg ratios by using the 
star particles identified to belong to the central spheroidal 
components. Those corresponding to the whole stellar mass 
and the total baryonic mass were also estimated and anal¬ 
ysed. These relations can be seen in Fig. [5] for z = 0 (upper 
panel), z = 1 (middle panel) and z = 2 (lower panel). As 
shown, there is a clear correlation for the disc and cen- 
tral spheroidal components in the three redshifts (se e also 
iDominguez-Tenreiro et al.lll998l : IPedrosa et al.l[20l3l . 

Let us first analyse the relations at z = 0. As can be 
seen in Fig. [21 there is a one-to-one relation between ji and 
mi, where i represents the gas or stellar discs components 
as expected in the case of angular momentum conserva¬ 
tion. As discussed before, the one-to-one correlation is the 
consequence of the fact that all the material in the system 
experienced the same external torques before it collapsed 
to form a bounded system, and that the collapse occurred 
with angular momentum conservation. Interestingly, for the 
stellar central spheroidal components a clear correlation is 
also detected but off-set from that determined by the disc 
components. This indicates that the material which ended 
up dominated by dispersion lost angular momentum as ex¬ 
pected, but it did so by a similar fraction, regardless of the 
fi nal stellar fract i on, m \. This trend wa s previously shown 
i niPedrosa et al.l (2014 figure 3) and in P^rosa fc Tisseral 


(120151 see also lDomfngue^Tenreiro et al.| (|l99fj ) for earlier 


results with lower numerical resolution and a different nu¬ 
merical code and subgrid physics). 

In Fig. 12] the same relations are estimated at z = 1 
and z = 2. As can be seen, both the stellar and gaseous 
discs are always consistent with the one-to-one relation. 
However, the most interesting result is again provided by 
the spheroidal components which determine a similar cor- 


Note that capital letters will be used for the total angular 
momentum and total mass and to denote the dynamical com¬ 
ponents: disc (D), central spheroid (S) or dark matter halo (H). 
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Fig. 1. Upper panels: Density map for a face-on and an edge-on view of a late type galaxy from the 230D simulation 
catalog. Lower panels: Disc and central spheroid components. 


relation as a function of redshift. Regardless of their forma¬ 
tion history, the disc and spheroidal components determine 
relations with the stellar and dark matter masses which are 
statistically unchanged up to z ^ 2. The gap between the 
angular momentum content of the galaxy respect to that 
of the dark matter does not show a systematic change with 
redshift. 

Hence, our simulated galaxies have disc and central 
spheroidal components consistent with the relation j; = 
A X to “, where A ~ 1 for the disc components and A ~ 
0.10 — 0.20 for the central spheroidal components. We fit 
a linear regression to the log-log relations finding slopes 
of a ^ 0.9 and a ^ 1.1 for the stellar discs and central 
spheroids, respectively at z = 0. 

These relations are found not to evolve significantly with 
redshift at least up to z ^ 2 (a ^ 0.90,0.87 and a ~ 
0.96, 0.86 for the stellar discs and spheroids at z = 1 and 
z = 2, respectively) Therefore, the spheroidal components, 
albeit dominated by dispersion, retain a similar fraction of 
the angular momentum imprinted by tidal torques onto the 
small perturbati ons during th e first stages of the evolution 
of the structure (|PeebleJl969ll . regardless of their assembly 
histories or mi. 

The slope of the relation ji-mi changes when the total 
stellar mass (i.e. total stellar discs and bulges) or the bary- 
onic mass are considered for the estimation of these ratios 
(blue and black lines, respectively), reaching values close to 
a 2 due to the combination of angular momentum con¬ 
tent of both disc and bulge components. This is in agree¬ 
ment with the results obtained by previous works where 
the angular momentum of baryons were studied without 
discerning betwe en disc and spheroidal components (e.g. 
ISales et ^11201211 . In this case, the gap is filled with the 


combination of contributions from both components and 
the slope of the mean relation gets steeper. We found indi¬ 
cations of a flattening in the relation for the total baryonic 
mass with redshift (a ~ 1.54,1.40 for z = 1 and z = 2, 
respectively ). However, a larger sample is needed to con¬ 
firm the level of evolution since the dispersion in our sample 
increases with increasing redshift. 

We have also analysed the angular momentum efficiency 
Tji for both components, defined as the ratio between jd (js) 
of the disc component (ce ntral spheroids) a nd jh of the 
dark matter halo, following ISales et ahl (l2niflh . In Fig.[3]we 
plot rji versus mi for the three analysed redshifts. As ex¬ 
pected, the disc components are found to have ^ 1, on 
average, reflecting the conservation of angular momentum 
during their formation. Conversely, the central spheroid has 
r] ^ 0.10 — 0.20, on average. In both cases, there is no 
clear evolution with redshift (although the dispersion in¬ 
creases slightly with increasing z). The different locations 
of r] for the two dynamical components reflect their different 
evolutionary histories which introduce dispersion but, nev¬ 
ertheless, take place in an organized fashion. The narrow 
range of mi values detected for the spheroids indicate that 
there are not tiny bulges or massive spheroids in our sim¬ 
ulations. Hence, our results should be verified using larger 
and higher-resolution simulations. 

3.1. The specific angular momenta of galaxies and their dark 
matter haloes 

In this section, we analyse in detail the j of the disc and cen¬ 
tral spheroids in relation to that of their dark matter haloes. 
In Fig.m we show j for the stellar disc and spheroidal com¬ 
ponents as a function of virial mass of their host galaxies. 
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Fig. 2. Relation ji — mi medians for the stellar (red dashed lines) and gaseous (green dashed lines) disc components, 
the stellar spheroidal components (pink dashed lines), the total stellar component (blue line) and the total baryonic 
component (black lines). The solid black line represents the relation in the case of angular momentum conservation. The 
stellar and gas disc components are consistent with angular momentum conservation while the spheroidal components 
determine a relation with a similar slope but off-set by a constant factor suggesting that spheroids have lost a similar 
fraction of angular momentum with respect to the total angular momentum of their virial haloes. The relation has been 
estimated for z = 0 (upper panel, in the right one we include the scatter data for the stellar disc component), and for 
z = 1 and z = 2 (bottom panel, left and right respectively), finding similar trends. 



-4 -3 -2 -1 0 

log m. 


Fig. 3. Angular momentum efficiency ryi as a function of 
m; ratio for disc (squares) and central spheroid (circles) 
components at z = 0 (red), z = 1 (green) and, z = 2 
(blue). 


They define parallel relations consistent with j oc 


This is in agreement with the expected theoretical relation 
derived under the hypothesis that the dark matter halo spin 
parameter is constant. For comparison purposes, we have 
also included the estimations for the gaseous disc compo¬ 
nents and the simulated dark matter haloes which agree 
well with the expected theoretical relation under the hy¬ 
pothesis of angular momentum conservation (solid line). 

Both, the spheroidal and disc components determine re¬ 
lations with the same slopes, which indicates that the stel¬ 
lar spheroids have formed from material of lower content 
of angular momentum (by a similar fraction), regardless of 
the virial mass of their host galaxy (or the relative frac¬ 
tion of stellar mass accumulated in the central region as 
shown in Fi g . [21). T his trend is similar to that reported by 
iGenel et al.l (1201511 . This agreement is greatly encouraging 
considering the differences in the codes, sub-grid physics 
and criteria adopted to classify the simulated galaxies. 

3.2. Evolution with redshift 

A question naturally arises about the possible evolution 
of j in the disc and central spheroidal components with 
redshift. This might imply a change in the slope or fraction 
of lost angular momentum, or both, as galaxies evolve. To 
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Fig. 4. Medians of the specific stellar angular momentum 
of the gaseous (green line) and stellar disc (red, dashed 
line) components and the central stellar spheroidal compo¬ 
nents (pink, dotted line) as a function of the virial mass of 
their host haloes. For comparison purposes, the relation es¬ 
timated for the dark matter haloes (black dotted line) and 
the expected theoretical relation (black solid line) have been 
also included. 

explore this, we estimated the j values for the stellar discs 
and central spheroids as a function of M^ir for z = 1 and 
z = 2. Recall that at each redshift, we apply the dynamical 
decomposition explained in Section 2.1. 

As shown in Fig. [Sj we found no significant evolution 
of these relations since z ^ 2. The gap between the j 
of rotational-dominated and dispersion-dominated compo¬ 
nents remains approximately constant, implying a lost of 
80 — 90% of the initial angular momentum content by the 
material which ends up in the spheroidal component, re¬ 
gardless of the virial mass of the halo or of the redshift. 

This is a remarkable trend which suggests again that, 
independently of the assembly history of galaxies, their disc 
and bulge components will contain similar amounts of an¬ 
gular momentum as they evolve. 

Note that in Fig. [5] we show the j for the stellar disc and 
bulge components, independently of the global morphology 
of the galaxies. Hence, some of the simulated galaxies are 
disc-dominated while others are spheroid-dominated sys¬ 
tems. So, their total j will be distributed filling the gap 
and determining a steeper mean relation as shown in in 
Fig. m Therefore, these trends of conservation of the j 
by rotational-dominated components and equally lost of j 
by dispersion-dominated ones will be erased as the gap is 
filled with the combination of systems. The separation of 
discs and central spheroids is crucial to get a deeper insight 
into the physical problem. Nevertheless, we can not explain 
the origin of the similar angular momentum efficiency for 
central spheroids yet. 

3.3. Comparison with observations 

Recently, I Fall fc Romanows^ (|2013ll reported results 
which support the existence of a relation between the j and 
stellar mass with the same slope f or both ellipticals and dis c 
dominated-galaxies, as found bv iFall fc Efstathi^ (|1980ll . 
The gap between them indicates a lost of angular momen- 


Fig. 5. Specific stellar angular momentum for the stellar 
disc (empty triangles) and central spheroid (empty circles) 
components as a function of the virial mass for z = 0 (red), 
z = 1 (green) and, z = 2 (blue). The gap between the two 
dynamical components does not evolve with redshift. For 
comparison purposes j for the the dark matter haloes have 
been also included (filled squares). Note that the angular 
momenta of the discs and central spheroids are given at 

?’HM- 


turn by the material which formed the ellipticals by a simi¬ 
lar fraction, regardless of stellar mass. Our simulations seem 
to reproduce naturally this trend as can be deducted from 
Fig. [5] and more clearly seen from Fig. E] where the o bserva- 
tional estimations from iFall fc Romanowskvl (1201311 are in¬ 
cluded. As shown in this figure, the stellar disc component 
determine a relation which agrees well with the one esti¬ 
mated for observed discs. In the case of central spheroids, 
they are also consistent with observations provided that 
stellar bulges behave similarly to elliptic al galaxies (as also 
suggested bv iFall fc Romanowskvl [201^ . The good agree¬ 
ment between the observations and simulations suggests 
that the latter are able to grasp well enough relevant as¬ 
pects of the physics at play. 

In order to improve the comparison of our results with 
observations, we estimated the total j defined by the whole 
stellar populations within the galaxy radiufl Previous fig¬ 
ures show results estimated at chm- The morphology of 
a simulated galaxy is defined by the ratio D/T between 
the stellar disc mass, D, and the total stellar mass of the 
galaxy, T = D + B. Simulated galaxies are classified accord¬ 
ing to the D/T in four subsamples from disc-dominated 
{D/T > 0.6) to spherodial-dominated {D/T < 0.2) ones. 

We should keep in mind that the comparison should be 
taken as indicative since differences between the observed 
and simulated parameters are always expected. For exam¬ 
ple, while observations use the gas information to trace the 
angular momentum content of the stellar discs under the 
hypothesis that they should be dynamical comparable, the 
simulated calculations make use of the information store 
directly in the stellar populations. As we show below, the 
gaseous discs conserve their angular momentum and they 
compare well to the relation determined by Sc galaxies, 

^ The galaxy radius is defined as the one which encloses ~ 80% 
of the baryonic mass. 
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where the dispersion-dominated component is small. This 
indicates that this type of galaxies is missing in the simu¬ 
lations. 

In Fig. m we show j as a function of galaxy stel¬ 
lar mass with different D/T ratios. A correlation can be 
seen in the sense that higher D/T ratios are related to 
higher contents of specific angular momentum, as expected. 
It also presents a good agreement w i th the observational 
trend found by iFall fc Romanows^ (l2013ll . The gap be¬ 
tween disc-dominated systems [D/T > 0.6) and spheroidal- 
dominated ones {D/T < 0.2) is filled with the correspond¬ 
ing intermediate morphological types. Our results support 
the claim that this relation opens the possibility to an al¬ 
ternative way of mapping galaxies’ morphologies. However, 
and more importantly, the agreement between the observa¬ 
tional and numerical results suggests that the hierarchical 
clustering signature is not erased in spheroidal galaxies as 
previously assumed. 

It is remarkable that the gap between both morpholog¬ 
ical types does not depend on the mass of the gala xy or 
the virial mass of the haloes. In iPedrosa et al.l (|2014ll a de¬ 
tailed study of the D/T fraction was done. The evolution 
of the morphology was quantified by analysing the D/T 
ratios as a function of the redshift for a group of 15 se¬ 
lected galaxies and the correlation of that parameter with 
their assembly history. As galaxies evolve, the D/T ratios 
of the progenitor systems can drastically change, reflecting 
their different formation histories (major and minor merg¬ 
ers, accretions, etc). The D/T fraction is a highly dynamic 
parameter, which indicates that the disc can be destroyed 
and re-buil t several times along its evolution, as shown in 
Figure 4 of iPedrosa et al.l (12014 ). 

During galaxy assembly, angular momentum is ex¬ 
changed thro ugh different phys ical processes. Previous 
works such as iBrook et al.l (1201 ll) have found that prefer¬ 
ential remotion of low angular momentum material at early 
redshift also contributes to the redistribution of the angu¬ 
lar momentum. They found that low angular momentum 
material is accreted and rapidly expelled when star forma¬ 
tion peaks at early times, with shallow potential wells that 
favour early gas ejection. Later accreted material has higher 
angular momentum and is accreted preferentially onto the 
disk, contributing to re-build it. Also, gas might flow into 
galaxies through the dark matter filaments, feeding the 
ga laxy with high an g ular momentum m aterial, as found 
bv iDekel et al.l ( 20091) . lOenel et all (l2015l) studied, using a 
more observational approach, the angular momentum con¬ 
tent of galaxies in the Illustris simulations. They found, in 
agreement with our results, that early-type galaxies reten¬ 
tion factor is approximately 30%. They also analised the 
evolution track of a group of selected galaxies, finding that 
different processes, such as mergers, could drastically redis¬ 
tribute the angular momentum content. 

Therefore, bulges are formed by important contribu¬ 
tions of low angular momentum material wh ich might have 
lost i t during the galaxy formation process ([Pedrosa et al.l 
120141) . 

4. Size relations 

As shown in previous sections, the disc components in our 
simulations were formed in agreement with the hypothesis 
of j conservation, while the central spheroids were formed 
from material with lower content of j, by a similar fraction. 



Fig. 6. Relation between j and the total stellar mass for 
the simulated stellar disc components (blue, filled squares) 
and the stellar spheroids (red, open circles). The observa - 
tional estimations reported by iFall fc Romanowskvl (|2013D 
are depicted for comparison: disc components (blue, open 
squares) and Elliptical galaxies (red filled circles). The the¬ 
oretical relation j oc is also included for comparison 

(black, continuous line). 


From these results, we would expect a clear scaling relation 
between the galaxy and halo size. 

The proportionality between the specific angular mo¬ 
mentum of the baryons and that of the halo assumed by 
the disc galaxy formation models would also imply the ex¬ 
istence of a linear relatio n between the g alactic and virial 
radius of the hosting halo (lKravtsovll2013l) . The galaxy size 
is expected t o scale as rgai oc Ar 2 oo • 

iKravtso^ (|2013l) analysed the relation between sizes of 
stellar systems of galaxies characterized by their half-mass 
radius and the virial radius of their dark matter haloes. 
The latter are derived by using the abundance matching 
ansatz. They found that the characteristic size of stellar 
and gas distributions in galaxies spanning several orders 
of magnitude in stellar mass, scales approximately linearly 
with a fraction of the the virial radius with a proportion¬ 
ality factor of ^ 0.015. However, this author did not con¬ 
sider the effects of baryons on the dark matter potential 
well, which are n aturally provided by our hydrodynamica l 
simulat ions ('e.g. iP edrosa et al.l l201fli : iTissera et al.l l2010l) . 
In fact. lDiemer et al.l (120131) claimed that due to the most 
common definition of the dark matter halo that considers a 
spherical overdensity of matter with respect to a reference 
density, such as the critical one or the mean matter density 
of the Universe, could lead to a spurious pseudo-evolution of 
halo mass due to redshift evolution of the reference density. 
They found that for all halo masses a significant fraction 
of the halo mass growth since z ~ 1 is due to this so-called 
pseudo-evolution rather than the actual physical accretion 
of matter. They obtained a corrected value for the size re¬ 
lation of 0.03 — 0.04. by taking the contraction effect into 
account. 

In Fig. |S] (upper panel) we show the stellar chm as a 
function of a fraction of r 2 oo for our simulated systems. 
Interestingly, we determined that the best factor relating 
both characteristic radii is ^ 0.04. In this figure, we show 
the relation for galaxies with different D/T ratios. It is clear 
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Fig. 7. Relation between j and the total stellar mass for simulated galaxies with different disc-to-bulge ratios: D/T > 
0.6 (red, filled triangles), 0.4 < D/T < 0.6 (green, filled triangles), 0.4 < D/T < 0.2 (blue, filled triangles) and 
D/T < 0.2 (pink, filled triangles) calculated at the gal axy radius. The relat i on for the gaseous discs are also included 
(cyan crosses).The observational estimations reported bv iFall fc R.omanowskvl (|201.lll are depicted for comparison: Sa-Sab 
(black triangles), Sb-Sbc (black squares), Sc-Sm (black diamonds), E (black filled circles), SO (black filled squares, sE 
(black circles). The theoretical relation j oc is also included (black, continuous line). 


that higher D/T values correspond to more extended struc¬ 
tures. And this is consistent with Fig.[7]where it can be seen 
that also higher D/T ratios correspond to higher values of 
angular momentum content. The size correlation is at place 
since z ~ 2 as expected since the angular momentum rela¬ 
tion is also present and unchanged globally since this time 
Fig. IHKlower panel). 


5. Conclusions 

We used hydrodynamical cosmological simulation to per¬ 
form a comparative analysis of the angular momentum con¬ 
tent in the disc and bulge components of the simulated 
galaxies. We focused our analysis on the the best resolved 
galaxies and used a dynamical criterion to identify the disc 
and bulge components at z = 0 and for their progenitor 
galaxies at z = 1 and z = 2. 

Our key findings are: 


We found that both discs and bulges determine a rela¬ 
tion in the ji — mi plane. Both the rotational-supported 
and the dispersion dominated components follow a re¬ 
lation with the same slope (j) oc mi). For the disc 
components, the correlation is the expected one-to- 
one relati on in the case of angular momentum con¬ 
servation (iMo et al.lll998ll . Remarkably, the dispersion- 
dominated systems determine a parallel relation but 
off-set by an approximately constant factor. Our sim¬ 
ulated bulges conserved ry ~ 10 — 20% of their origi¬ 
nal angular momentu m, regardless of the D/T ratio of 
the galaxy (see als o Domfnguez-Tenreiro et al.l 119981: 
iPedrosa et al.l[2014 ). This is in agreement with observa¬ 
tional findings ( Fall &: Romanowskvll201^ and referen- 
cies there in). 


— A similar analysis was carried out for the progenitor 
systems up to z ~ 2, finding that the central spheroids 
and discs determine similar relations. Hence, regardless 
of the formation history of these galaxies, their disc and 
spheroidal components have angular momentum con¬ 
tents which determine relations with the stellar and 
dark matter masses with similar slope up to z 2. In 
the case of the spheroidal components a similar off-set 
is detected showing an angular momentum efhciency of 
^ 0.10 — 0.20. These values vary slightly with redshift. 
A larger sample is required to assess robustly the level 
of evolution. 

— Regarding the specific angular momentum, we found 

that the spheroidal and disc components determine re¬ 
lations with the same slopes, j oc This indicates 

that the stellar spheroids have been formed by material 
of low angular momentum content or that have lost it by 
similar fraction during the assembly process, regardless 
of the virial mass of their host galaxy. We have analysed 
this relation at different redshifts, finding no evolution. 
Remarkably, this trend is suggesting that independently 
of virial masses, the disc and bulge components will con¬ 
serve similar amounts of angular momentum as they 
evolve. 

— In good agreement wi th the trend found by 
iFall fc Romanowskvl(|2013ll . we found a clear correlation 
between the morphological type of a galaxy and its total 
specific angular momentum content: higher D/T ratios 
are related with higher contents of specific angular mo¬ 
mentum, as expected. The gap between disc-dominated 
systems {D/T > 0.6) and spheroidal-dominated ones 
{D/T < 0.2) is filled with the corresponding intermedi¬ 
ate morphological types. This supports the claim that 
this relation opens the possibility to an alternative way 
of mapping the morphology of galaxies. And this also 
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Fig. 8. Size relations for the simulated galaxies and their 
host dark matter haloes: thm = 0.04 x Ar 2 oo. Upper panel: 
Simulated galaxies are colour-coded according to the D/T 
ratio: D/T > 0.6 (red, filled squares), 0.4 < D/T < 0.6 
(green, filled diamonds), 0.4 < D/T < 0.2 (blue, crosses) 
and D/T < 0.2 (pink, filled circles). Lower panel: The 
colours correspond to the different redshifts, z = 0 (red 
circles), z = 1 (green squares), z = 2 (blue triangles) 


implies that the hierarchical clustering signature is still 
present in spheroidal galaxies. It is remarkable that the 
gap between both morphological types does not depend 
on the mass of the galaxy or the virial mass of the 
haloes. 

— In agreement with lKravtso\l (l20I3ll , we found that the 
characteristic size of stellar distributions in the simu¬ 
lated galaxies scales approximately linearly with the 
virial radius in the log-log plane. The p roportionality 
facto r for our simulated systems is ^ 0.04 (jDiemer et al.l 
1201,111 . This relation is found not to evolve with time: 
galaxies move along the relation as one moves to higher 
redshift. 
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